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I. INTRODUCTION 



Recently there has been considerable progress both in the visualization (Bednar et al., 1998) and micromanipulation 
(Cui and Bustamante, 2000) of chromatin fibers. These results constitute an important step towards the understanding 
of DNA "folding", i.e., the problem of how plant and animal genomes organize themselves into volumes whose linear 
dimensions are many orders of magnitude smaller than their contour lengths. For instance, human DNA is billions of 
base pairs (bp) long (about a meter), this length of highly-charged (about one fundamental charge per two Angstroms) 
and hard-to-bend (persistence length of 50nm) linear polymer must be condensed into chromosomes that fit into cell 
nuclei whose characteristic size is a micron. 

An important part of the condensation process is the complexation of DNA with oppositely-charged globular protein 
(histone) aggregates that have the shape of squat cylinders. These aggregates are octameric complexes consisting of 
pairs of the four core histones H2A, H2B, H3, and H4. A DNA stretch of 147 bp's is wrapped in a 1- and- 3/4 
left-handed superhelical turn around the histone octamer and is connected via a stretch of " linker" DNA to the next 
such protein spool. Each protein aggregate together with its wrapped DNA comprises a nucleosome core particle (cf. 
Fig. 1) with a radius of about 5nm and a height of about 6nm; with its linker DNA it is the fundamental chromatin 
repeating unit. It carries a large negative electrostatic charge (Khrapunov et al., 1997; Raspaud et al., 1999). Whereas 
the structure of the core particle has been resolved up to high atomic resolution (Luger et al., 1997), there is still 
considerable controversy about the nature of the higher-order structures to which they give rise. When stretched, the 
string of DNA/histone complexes has the appearance of "beads-on-a-string" . This basic structure can be seen clearly 
when chromatin is exposed to very low salt concentrations, and is known as the 10-nm fiber (Thoma et al., 1979), 
since the diameter of the core particle is lOnm. With increasing salt concentration, i.e., heading towards physiological 
conditions (lOOmM), this fiber appears to thicken, attaining a diameter of BOnm (Widom, 1986). The absence of the 
extra "linker histones" (HI or H5) leads to more open structures (Thoma et al., 1979) so it is surmised that the linker 
histones act near the entry-exit point of the DNA (cf. Fig. 1); they carry an overall positive charge and seem to 
bind the two strands together leading to a stem formation (Bednar et al., 1998). Increasing the salt-concentration is 
expected to decrease the entry-exit angle of the stem as it reduces the electrostatic repulsion between the two strands. 

Longstanding controversy (van Holde, 1989; Widom, 1989; van Holde and Zlatanova, 1995, 1996) surrounds the 
structure of this 30-nm fiber, for which there are mainly two competing classes of models: the solenoid models (Finch 
and Klug, 1976; Thoma et al., 1979; Widom and Klug, 1985); and the zig-zag or crossed-linker models (Woodcock et 
al., 1993; Horowitz et al., 1994; Leuba et al., 1994; Bednar et al., 1998). In the solenoid model (Fig. 2a) it is assumed 
that the chain of nucleosomes forms a helical structure with the axis of the core particles being perpendicular to the 
solenoid axis (the axis of an octamer corresponds to the axis of the superhelical path of the DNA that wraps around 
it). The DNA entry-exit side faces inward towards the axis of the solenoid. The linker DNA (shown as a dashed curve 
at the top of Fig. 2a) is required to be bent in order to connect neighboring nucleosomes in the solenoid. The other 
class of models posits straight linkers that connect nucleosomes located on opposite sides of the fiber. This results in 
a three-dimensional zig-zag-like pattern of the linker (Fig. 2b). 

Images obtained by electron cryomicroscopy should in principle be able to distinguish between the structural features 
proposed by the different models mentioned above (Bednar et al., 1998). The micrographs show a zig-zag motif at 
lower salt concentrations and they indicate that the chromatin fiber becomes more and more compact when the ionic 
strength is raised towards the physiological value. However, for these denser fibers it is still not possible to detect the 
exact linker geometry^. 

An important experimental achievement was the stretching of a single chromatin fiber via micromanipulation (Cui 
and Bustamante, 2000). The "force-extension" measurements show a rich behavior of the mechanical properties 
as a function of the ionic strength. At low ionic strength (5raM NaCl) the force-extension curves are reversible 
as long as the tension does not exceed 20pN. For higher tension levels (> 20pA^) there are irreversible changes 
that lead to an increase of the fiber length, probably due to the loss of linker histones and/or histone octamers. 



^Experiments on dinucleosomes (two nucleosomes connected by one linker) have been performed to check if the nucleosomes 
"collapse" upon an increase in ionic strength. A collapse would only occur if the linker bends, and an observation of this 
phenomenon would support the solenoid model. The experiments by Yao et al. (Yao et al., 1990) as well as more recent 
experiments by Butler and Thomas (Butler and Thomas, 1998) indeed reported a bending of the linkers but do not agree 
with experiments by Bednar et al. (Bednar et al., 1995) and by others that did not find any evidence for a collapse. Critical 
discussions of these and other experiments on dinucleosomes are available (van Holde and Zlatonova, 1996; Widom, 1998). 
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At high ionic strength (40mM and 150mM NaCl) a 5piV-plateau in the force extension curve was found^. The 
authors interpreted their results as indicating a coexistence between "swollen" and "condensed" parts of the fiber. 
In order to reproduce the force-extension curves, Katritch et al. performed Monte Carlo simulations (Katritch et al., 
2000) based on a geometrical "two-angle" model introduced by Woodcock et al. for the 30-nm fiber (Woodcock et 
al., 1993) combined with the worm-like chain ("WLC") Hamiltonian for the linkers. The WLC is widely used for 
predicting the mechanical properties of naked DNA. The low-salt behavior could be reproduced for several sets of 
angles and bond lengths of the model, demonstrating both that it is a reasonable model and that it is not possible 
to deduce a unique structure from the measured response of the fiber under stretching. Katritch et al. found that an 
internucleosomal attraction of roughly 3kT might explain the experimentally observed plateau in the force-extension 
profile. The biological importance of these results lies in the fact that significant changes can be achieved in the degree 
of chromatin condensation with only modest levels of external stress. The fact that chromatin at physiological salt 
concentrations apparently can exist in two alternative forms that interconvert under low levels of stress is particularly 
interesting. 

The success of the model motivated the present study to provide an analytical framework for understanding the 
geometrical and mechanical properties of the 30-nm fiber based on the two-angle model. Our first main result is 
the derivation of a general structural phase diagram of the chromatin fiber as a function of the two basic angles 9 
and (j) determined by the nucleosome properties and the linker length b (see below). The various solenoidal, zig-zag 
and crossed- linker structures - all of which are assumed to have straight linkers - appear as " points" in this phase 
diagram (see Fig. 4) . We find that, within the two-angle model, the position of chromatin fibers at physiological 
conditions (the "native" fibers) in the phase diagram is surprisingly close to the point in the diagram with the highest 
density and the maximal accessibility, consistent with excluded-volume restrictions. Changes in bond angles induced 
by physico-chemical changes in the environment lead to predictable changes of the fiber away from this optimal point 
towards more open structures. 

Our second main result is that we can obtain approximate analytical results for the bending stiffness of the two-angle 
model - and hence for the persistence length - and for the force-extension curve. We find (?) that the persistence 
length of the fiber should be comparable or less than that of naked DNA, for a wide range of 6- and ^-values and (ii) 
that the stretching modulus should be so low that there is no longer a pronounced difference between " soft" entropic 
elasticity (for low forces) and "hard" entropic elasticity (for high forces), in marked contrast with the case of naked 
DNA. Using the estimated values of 0, (j) ^^'^ ^ the predicted force-extension curves (with no fitting parameter) are in 
good agreement with the data found for the stretching of chromatin fibers (Cui and Bustamante, 2000). 

The implication of our results is that a swollen 30-nm fiber should be very soft in terms of its elastic properties, over 
a wide range of values of the angle parameter 6 and 0, a very reasonable "design feature" in terms of its biological 
role. This swollen state competes with a more rigid condensed state that appears, as a function of bond angle 6, 
when we allow for (weak) attractive forces between nucleosomes. The physical properties of the condensed state are 
beyond the scope of the current paper, but the condensed fiber is expected to be significantly stiffer that the swollen 
fiber. In general, our results appear to indicate that the "engineering design" of the 30-nm fiber combines high 
compaction levels with high structural accessibility and flexibility. Independent of the question whether the swollen 
or the condensed state is realized, modest changes in the control parameter tt — 9 (the nucleosome entry-exit angle) 
produce large structural changes. 

The paper is organized as follows. In the Section 2 we derive the geometrical properties of the two-angle model 
and present the general diagram of states. In Section 3 we apply our results of the two-angle model to interpret 
the structure of the 30-nm chromatin fiber in terms of simple optimization principles. Section 4 derives the elastic 
properties of the two-angle model and gives the bending stiffness and the force-extension relation. In the concluding 
section we summarize our results and discuss alternative models. 

II. THE TWO-ANGLE MODEL: FOLDED STRUCTURES 



■^Marko and Siggia (1997) had in fact proposed an elastic model which predicted a coexistence regime in the force-extension 
curve, with nucleosomes "evaporating" from the fiber at higher force levels of the order of 2pN which would lead to extensive 
irreversibility in the force-extension curve. Although irreversibility is encountered at high force levels, as mentioned, the 5pN 
plateau is reversible indicating that there was no nucleosomal loss. 
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A. General relationships 



To address the folding problem of DNA at the level of the 30-nm fiber we need a mathematical description for the 
different possible folding pathways. At the simplest level, it is assumed that the geometric structure of the 30-nm fiber 
can be obtained from the intrinsic, single- nucleosome structure. The specific roles of linker elastic energy, nucleosome- 
nucleosome interaction, preferred binding sites, HI involvement, etc. will be treated afterwards as "corrections" to 
this basic model. To see how single-nucleosome properties can control the fiber geometry, consider the fact that DNA 
is wrapped a non-integral number of turns around the nucleosome, e.g., l-and-3/4 times (147 bp's) in the case of 
no HI. This implies that the incoming and outgoing linker chains make an angle 9 with respect to each other — 
the entry-exit angle tt — is nonzero. In the presence of the histone HI (or H5) the in- and outcoming linker are 
in close contact forming a stem before they diverge (Bednar et al., 1998). While the precise value of the resulting 
exit-angle depends on salt concentration, degree of acetylation of the histones, etc., we may nevertheless assume 9 
to be a quantity that is determined purely at the single-nucleosome level. Next, we define the rotational (dihedral) 
angle (p between the axis of neighboring histone octamers along the necklace (see Fig. 3). Because nucleosomes are 
rotationally positioned along the DNA, i.e., adsorption of DNA always begins with the minor groove turned in towards 
the first histone binding site, the angle is a periodic function of the linker length b, with the lObp repeat length of 
the helical twist of DNA as the period. There is experimental evidence that the linker length shows a preferential 
quantization involving a set of values that are related by integral multiples of this helical twist (Widom, 1992), i.e., 
there is a preferred value of 0. (Note that the "linker length" b is strictly speaking defined here as the distance 
between two neighboring nucleosomes, cf. Fig. 3.) 

If we treat the pair of angles {9, 0), together with the linker length 6, as given physical properties (even though in 
vivo they are likely under biochemical control), then the geometrical structure of the necklace is determined entirely 
by 9, (j) and b. The model only describes linker geometry and does not account for excluded volume effects and other 
forms of nucleosome-nucleosome interaction; it assumes that the core particles are pointlike (a = 0) and that they are 
located at the joints of the linkers. The model also assumes that the linkers are straight. It is under dispute whether 
this last condition holds for the 30-nm fiber at higher salt concentrations, and we will return to this issue later. The 
(9, (/))-model is similar to the freely rotating chain model encountered in polymer physics literature (see, for instance, 
Doi and Edwards, 1986). The main difference is that in the present case there is no free rotation around the linker 
and so torsion is transmitted (see also Plewa and Witten, 2000). 

As shown in Appendix A it is now possible to construct a spiral of radius R and pitch angle 7 such that the 
nucleosomes - but not necessarily the linker chain - are located on this spiral. The nucleosomes are placed along 
the spiral in such a way that successive nucleosomes have a fixed (Euclidean) distance b from one another. From 
straightforward geometrical considerations we can derive analytical expressions that relate pitch angle 7 and radius 
R of the solenoid to the pair of angles 9, cj) and linker length b. Specifically the linker length b can be expressed as 
a function of 7, R and sq (defined as the vertical distance between successive "nucleosomes" along the helical axis), 
b — fe(7,i?. So), as given by Eq. p% The corresponding relationships for the angles 9 and (j), 9 = 9{'~f,R,so) and 
(/) = (/) (7, R, So), are given by Eq. pSfand Eq. ^ Using these relations, we can construct a catalog of structures. 

B. Planar structures 

If either one of the angles 9 oi (j) assumes the value or tt, then the resulting structure is planar^ and calculation 
of the associated geometrical properties is straightforward. Let us start with the case (j) ~ 0. If we also have 9 — 
the fiber forms a straight line (see structure "1" in Fig. 4). For small non- vanishing 9 the structure forms a circle of 
radius R b/9. For the special case 9 — 27r/n, with n an integer, the ring contains n monomers before it repeats 
itself and we obtain a regular polygon (see "2"). The special case is 9 — it/ 2 corresponds to the square ("3"). With 
increasing 9 the radius of the circle shrinks and approaches asymptotically the value 6/2. For 9 = tt [n — 1) jn with n 
being an odd integer one encounters a series of closed star-like polygons with n tips. In particular, n = 3 corresponds 
to the equilateral triangle ("4"), n = 5 to the regular pentagram ("5"), etc. 

Next we consider the case (j) — tt and 9 arbitrary. This case corresponds to 2D zig-zag-like structures, as shown by 
"6" and "7" at the top of Fig. 4. The length of a fiber consisting of N monomers is given by 

L ^bcos{9/2)N (1) 

and the diameter is given by _D = 6 sin (9/2). Note that the length of the fiber increases with decreasing 9. 
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To complete our discussion of planar structures we mention the remaining cases: = with an arbitrary value of 
4> leads to the straight line mentioned earlier ("1"); 9 = tt and arbitrary cj) corresponds to linkers that go back and 
forth between two positions (" 8" ) . 



C. Three-dimensional fibers 



(i) Solenoids: For small angles, 9 <^ 1 and ((> <ti 1, we find structures that resemble solenoids where the linkers 
themselves follow closely a helical path (see "9" in Fig. 4). For these structures one has aso/R <C 1 — where a = cot 7, 
with 7 the pitch angle. To lowest order in usq/R we find b ~ soVl + (cf. Eq. |2|), 9 ~ a^b/ (i? (l + a^)) (cf. 
Eq. |3^), and (j) ~ 9/a (cf. Eq. ^). From this we can infer several geometrical properties of the fiber as a function 
of b, 9 and (p summarized in Tabic I. R denotes the radius of the fiber, L is the length of a fiber consisting of + 1 
monomers, A is its line density N/L and p is the 3D density given hy p = X/ttR^, assuming a hexagonal array. 

Other geometrical information can be obtained easily. For instance, the vertical distance d between two loops 
follows from L in Table I by setting N = 2t:/9 (the number of monomers per turn): 



1^^T9^ 



(2) 



Furthermore, the pitch angle 7 is given by 



cot 7 ~ — (3) 



7 decreases monotonically as the ratio of the angles, 9/(f), increases. For 4> 9 one finds 7 ~ (j)/9. In this regime one 
has a very dense spiral with d ^ R. In the opposite limit (j)^ 9 the pitch angle is very large, namely 7~7r/2 — 0/0 
and the solenoid has a very open structure with d '3> R^^ 

(a) Fibers with crossed linkers: Consider structures where (f) is still small but where the entry-exit angle 9 is large, 
i.e. TT — 9^TT. We discussed in the previous section that for (/> = one encounters star-shape polygons that are closed 
for 9 = 7r(n — l)/n with n odd. For no?i- vanishing the star-shaped polygons open up in an accordion- like 

manner. This leads to a three-dimensional fiber with crossed linkers - see " 10" . It follows from Eqs. ^ and ^ that 
s§ ~ <jp (4i?^ — fo^) /4 for ^ 1. Using this result as well as Eq. ^ _R, L, A and p can be expressed as a function of 
b, 9 and 0, cf. Table I. 

Assume now that 6'„ = tt (n — 1) /n so that the projection of the fiber is a closed polygon (this is only strictly true 
for (p = but it is still a good approximation for <C 1). We can calculate for this case the spacial distance d between 
nucleosome i and i -\- n: 

d^^^ot{9/2)^^ 1 + —^] (4) 



2 ' ' ' 4 V 12" , 

{Hi) Twisted zig-zag structures: Finally, we discuss structures with a rotational angle close to tt, say (/) = t: ~ 5 
with (5^1. For S = we recover the 2D zig-zag structure discussed earlier ("6" and "7"). Small non- vanishing 
values of 6 lead to twisted zig-zag structures - see "11". In this case monomer i -I- 1 is located nearly opposite to the 
ith monomer, but slightly twisted by an angle S. Monomer z -I- 2 is then on the same side as monomer i but slightly 
twisted by an angle 26 and so on. The geometrical properties are given in Table I. For = tt, i.e., 6 — 0, we recover 
the result for the planar zig-zag structure. 

The fiber is contained within a cylinder of radius R, given in Table I. The monomers ("histones") are located at 
the surface, with the linker passing back and forth (approximately) through the middle axis of the cylinder. The 
monomers n, n ± 2, n ± 4... and the monomers n ± 1, n ± 3, ... form a double helix that winds around the cylinder. 
Within each of the two spirals the monomers are not directly linked together, even though monomer i and i + 2 can 
come quite close in space for large values of 9. The pitch angle of the two spirals follows from the positions of monomer 
1 and 3; Pi = (i?,0,0) and P3 ~ {R,-2R5,2bcos{9 /2)), cf. Eq. |3^. Thus 7 = -tt + A7 with A7 ~ (5 tan (61/2) /2. 



^We note that such an open structure could in principle collapse into a very dense fiber like the solenoidal model proposed by 
Klug (cf. Fig. 2a) if we would allow the linkers to bend. As mentioned already before it is still a matter of controversy if such 
linker bending takes place in chromatin. We will stick in this study to the assumption of straight linkers. 

*We mention that in the limit — > we recover the planar circle with radius Rc^h/O, ci. Table I. 
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D. Structure diagram and excluded volume restriction 



We now consider the full range of states in the {9, 0)-space as shown in Fig. 4. Both angles 6 and </> can each 
vary over the range to tt. At the edges of the diagram where one of the angles assumes an extremal value, the 
configurations are always planar. On the line = we find circles and star-type polygons (that are closed for specific 
values of 0). The planar zig-zag-structures are located on the line (j) = it; for = we find a straight configuration 
and for = tt a "dimer" structure. If we move from the line = towards larger values of 4> the circles and star-like 
polygons stretch out into the direction perpendicular to their plane, forming a solenoid and a fiber with crossed linkers, 
respectively. On the other hand, if we start at the top of the diagram (0 = tt) and decrease the value of 4> the planar 
zig-zag structure extends into the third dimension by becoming twisted. If we start with a structure with entry-exit 
angle = and increase the value of this angle, then the structure folds either into a solenoid with large pitch angle 
for small 0- value or into a twisted zig-zag for large values of </>. Finally, starting out at the dimer configuration, 9 = tt 
leads to an unfolding of the structure into a fiber with crossed linkers (small 0-values) or a twisted zig-zag (large 
(/)- values) . 

If we take into account the excluded volume of the core particles, then certain areas in our phase diagram are 
forbidden - reminiscent of the familiar " Ramachandran plots" used in the study of protein folding (Stryer, 1995). For 
simplicity we assume in the following that the core particles are spherical with a radius a and that their centers are 
located at the joints of two linkers, cf. Fig. 3. There are two different types of interactions. One is between monomers 
at position i and i ± 2 (short range interaction) , and leads to the requirement that the entry-angle must be sufficiently 
small: 

9 < 2 arccos (a/6) ~ tt ^ , a -4^ b (5) 

This condition excludes a vertical strip at the right side of the diagram, as indicated in Fig. 4 by a dashed line. 

There is also a long-range excluded volume interaction that comes into play when the angle 4> is too small. This 
is apparent for the case = where we find planar structures that run into themselves. Starting with a circular 
structure we have to increase (p above some critical value so that the pitch angle of the resulting solenoid is large 
enough so that neighboring loops do not interact. This leads to the requirement d > 2a with d given by Eq. |^ (using 
(j> < 9), i.e., 

0>-?^' (6) 
For the large 9-case (fibers with crossed linkers) we find from Eq. ^ the condition 

The two conditions, Eqs. |6| and ^ shown schematically as a dotted curve in Fig. 4, lead to a forbidden strip in the 
structure diagram for small values of (j). 

Figure 4 does not show the interesting "fine-structure" of the boundary of the forbidden strip that is due to 
commensurate-incommensurate effects. We already noted that there are special 6'-values for which the projection of 
the linkers forms a regular polygonal star (0„ = tt (n — 1) /n) or a regular polygon (61^ = 27r/n) (for small values of 
(j)). In these cases the nucleosomes i and i + n "sit" on top of each other. On the other hand, for other values of 9, 
monomers of neighboring loops will be displaced with respect to each other. In this case monomers of one loop might 
be able to fill in gaps of neighboring loops so that the minimal allowed value of cj) is smaller than estimated above. We 
have not explored the interesting mathematical problem of the exact boundary line since this is likely to be sensitive 
to the exact nucleosome shape. The dotted line in Fig. 4 only represents the upper envelope of the actual curve. 

Our discussion of the two-angle model was based of the assumption of a perfectly homogeneous fiber where 6, 9 
and (f> are constant throughout the fiber. The effect of randomness in these values on the fiber geometry is discussed 
in Appendix B. 



III. CHROMATIN AND THE TWO-ANGLE MODEL: OPTIMIZATION OF DESIGN? 

Where in the structure diagram is actual chromatin located? The classical solenoid model of Finch and Klug (Finch 
and Klug, 1976) is found in the small 9, small (p section of the diagram (although in their case the linker is bent). 
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Various structures were displayed by Woodcock et al. in their Fig. 2 (Woodcock et al., 1993), namely fibers with 

= 150° and many different values of (j>, corresponding to a vertical trajectory on the right-hand side of Fig. 4. Three 
different configurations with a fixed value of 4> and different values of are displayed in Fig. 3(c) in another paper by 
these authors (Bednar et al., 1998). 

Our structure diagram accommodates all of these structures and, by itself, does not favor one over another. However, 
our diagram plus the formulae given above are useful if we invoke the following two criteria to optimize the structure 
of the 30-nm fiber: 

(i) maximum compaction 

(ii) maximum accessibility 

The first criterion is obvious: inactive chromatin should be packed as dense as possible because of the very large ratio 
of DNA length to nucleus size. By the second criterion we mean that a local accessibility mechanism is required for 
gene transcription. 

In order to attain maximum compaction we need structures that lead to high bulk densities p (we assume that 
the 30-nm fibers are packed in parallel forming a hexagonal lattice). A comparison of the 3D densities of the three 
different structures given in Table I shows that fibers with internal linkers have highest densities p, namely 

^~ TT(f>{TT~0)b3 ^ ^ 

In particular, the highest density is achieved for the largest possible value of 9 and the smallest possible value of 
(j) that is still in accordance with the excluded volume condition. This set of angles is located at the point where 
the dotted curve and the dashed line in Fig. 4 cross each other. Apparently this also represents the only region in 
the phase diagram where excluded volume effects are operative on a short-range and a long-range scale at the same 
time, i.e., nucleosome i is in close contact with nucleosome i — 2 and i + 2 as well as with nuclesomes father apart 
along the contour length of the necklace. This unique set of angles is given by 6'max ~ 2 arccos (a/5), cf. Eq. ||, and 
0min ~ (8/7r)(a/6), cf. Eq. 

In order to achieve maximum accessibility we look for structures that, for a given entry-exit angle tt — of a highly 
compacted structure, achieve the maximum reduction in nucleosome line density A for a given small change A9 of 
the angle 0. In other words, we look for a maximum of dA/d^ which we call the "accessibility". Interestingly, the 
accessibility is maximized at the same unique pair of angles (0max,</'min)- This can be seen from its angle dependence 
for fibers with crossed linkers 

! (9) 

d0 (j)h{-K-ef 

We note that this change in A with 6 is achieved by changing the number of monomers per vertical repeat length d. 
The length d itself is only weakly dependent on n according to Eq. ^. 

Before we compare our theoretical formulas with experimental results we mention that for fibers with crossed 
linkers there might be another excluded volume interaction, namely between linkers. For these structures the linker 
connecting the monomers i and i + 1 comes closest to the linkers between monomer i + 2 and i -I- 3 and the one between 

1 — 1 and i — 2, as can be seen for the n = 5 case, cf. "5" and "10" in Fig. 4. The linkers cross close to the middle of 
the fiber where the distance between their axes is given by 2 {4)b/2) cot [0/2) ~ 05 (tt - 0) /2, cf. Eq. |. This distance 
is minimized at (0max,</'min) and has to be larger than the thickness t of the fiber: 

We compare now the above given formulas with experimental results. For chicken erythrocyte chromatin one has 
roughly h ~ 20nm (center-to-center distance of nucleosomes, van Holde and Zlatanova, 1996). Together with a « 5nm 
this leads to 6'inax ~ 151°, (j>min ~ 36° and A w 6.9 nucleosomes per llnm (cf. Eqs. ||, and ||). Furthermore, the 
condition on the linker thickness is given by t < 3.2nm and is fulfilled since t = 2nm for DNA. The theoretically 
derived values can now be compared with the ones reported by Bednar et al. for chicken erythrocyte chromatin fibers 
(Bednar et al., 1998). From their table Iwe find that for an ionic strength of 80mM (which is close to the physiological 
value) « 145° and A = 5.9 nucleosomes per llnm. Furthermore, electron cryotomography-constructed stereo pair 
images of an oligonucleosome (cf. Fig. 3(b) in Bednar et al., 1998) indicate that the chromatin fiber might indeed 
have the structure of a fiber with crossed linkers, with ti « 5; this would correspond to 6* = tt (n — 1) /n w 144°. 



7 



Information concerning the preferred value for cj) may be obtainable from the measured statistical distribution of 
the nucleosome repeat lengths. This distribution shows statistically preferred linker lengths equal to lOfc + Ibp's with 
k a positive integer (Widom, 1992), which, in turn, indicates that the rotation angle (j) corresponds to a change in 
helical pitch associated with Ibp, i.e. 360°/10 = 36°. This value coincides with 0min, the value that we estimated for 
maximum compaction.^ 

The second feature, the local accessibility, can be monitored in vitro by changing the salt concentration. Bednar 
et al. report, for example, that 9 decreases with decreasing ionic strength, namely 6 ^ 145° at 80mM, 9 ~ 135° 
at 15mM and 9 ~ 95° at 5mM (Bednar et al., 1998). In the biochemical context the change of 9 is accomplished 
by other mechanisms, especially by the depletion of linker histones and the acetylation of core histone tails (cf., for 
instance, van Holde and Zlatanova, 1996), both of which are operative in transcriptionally active regions of chromatin. 
These mechanisms lead effectively to a decrease of 9. 

As pointed out below Eq. ^, the decrease of 9 is accompanied by a decrease of the line-density A = n/d oi 
nucleosomes at an essentially fixed value of d. In other words, the number of vertices of the projected polygon 
decreases significantly with decreasing 9 because 0„ = 7r(l— 1/n). In that respect the effect of reducing 9 below 
the optimal packing value might be best viewed as an "untwisting" of the 30- nm fiber. Using the experimentally 
determined values of 9 we find from Table I that the density (the number of nucleosomes per llnm) is given by 
A « 6.8 for 9 « 145°, A « 4.5 for 9 « 135° and A « 2.3 for 9 « 95°, slightly higher than the experimental values 
A « 6.0, A w 3.2 and A « 1.5 (Bednar et al., 1998). Furthermore, the number of polygonal vertices n = tt/ (tt ~ 9) 
decreases as follows: n « 5.1 for 9 « 145°, n « 4.0 for 9 « 135° and n « 2.1 for 9 « 95°, consistent with the stereo 
pair images by Bednar et al., suggesting n « 5 at an ionic strength of 80mM and n « 3 at 5mM (cf. Figs. 3(a) and 
(b) in Bednar et al., 1998). 

We close this section with a cautionary remark. The 3D density and the line density of the fiber can not only be 
changed by changing 9 or 4> but also by changing the linker length (in multiples of lObp's). A variation in b changes 
the location of the point (6'max, 0min) in the diagram of geometrical states, and thus the values of the maximum 3D 
and line densities that can be achieved, namely 

(11) 



TT^min (tT ^ 6'max) a?h 

and 



TT 6 , , 

-- (12) 



6(/'min (tt - 0max) 4 a 

This shows that fibers with smaller values of b can achieve higher 3D densities but have a smaller maximal line 
density (and accessibility dX/d9 esc b^). From this one might infer that active cells should have larger nucleosome 
repeat lengths in order to maximize the accessibility to their genetic material. An overview on nucleosome repeat 
lengths in different organisms and tissues is given in table 7-1 of van Holde's book (van Holde, 1989). The data shown 
there do not follow this rule, unfortunately. In fact, very active cells like yeast cells and neuronal cells have in general 
short nucleosome repeat lengths while inactive ones like sperm cells have large ones. This shows that the optimization 
principle of high density has to be used with caution. 



IV. STRETCHING AND COMPRESSION OF TWO-ANGLE FIBERS 



A. Introduction 



The (0, 0, b) model developed so far is purely geometrical. Could it be useful as well for predicting physical properties 
of the 30-nm fiber? The response of the 30-nm fiber to elastic stress will be the focus of this section. The elastic stress 
can either be of external or of internal origin. External stresses are exerted on the chromatin during the cell cycle 
when the mitotic spindle separates chromosome pairs. The 3Q-nm fiber should be both highly flexible and extensible 



^The statistical uncertainty around the expectation values for the nucleosome repeat length is sufficiently large to make our 
estimate for less reliable. 
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to survive these stresses. The in vitro experiments by Cui and Bustamante demonstrated that the 30-nm fiber is 
indeed very "soft" (Cui and Bustamante, 2000). 

The 30-nm fiber is also exposed to internal stresses. Attractive or repulsive forces between the nucleosomes will 
deform the linkers connecting the nucleosomes. For instance, electrostatic interactions, either repulsive (due to the 
net charge of the nucleosome core particles) or attractive (bridging via the lysine- rich core histone tails; Luger et al., 
1997) could lead to considerable structural adjustments of the (6*, (f>, b) model. 

In this section we will derive an analytical description of the force extension curve of the {6, 0, b) model in order to 
predict the elastic properties of the different structures obtained in the previous section. Using the particular values 
of and (j) that are observed experimentally (Bednar et al., 1998; Widom, 1992), we can reproduce rather well the 
measured force-extension curve of Cui and Bustamente and the numerical results of Katritch et al. that were based 
on a variant of the (6*, (j), b) model (see below). 

Before considering the elastic properties of the (6*, 0, b) model, it is helpful to briefly recall some results concerning 
the large-scale elasticity of the DNA itself (Cluzel et al., 1996; Marko, 1998). The measured force-extension curve 
of naked DNA breaks up into two highly distinct regimes: the "entropic" and "enthalpic" elastic regimes. For very 
low tension / (< \pN), the restoring force is provided by "entropic elasticity" (de Gennes, 1979). In the absence of 
any force applied to its ends, the DNA's rms end-to-end distance (chain length, L) is small compared to its contour 
length (Lo) and the chain enjoys a large degree of conformational disorder. Stretching DNA reduces its entropy and 
increases the free energy. The corresponding force / increases linearly with the extension L: 

/-f^f, L«L, (13) 

^DNA 



The length ^dna is known as the " thermal persistence length" of DNA and is of the order bQnm (Hagerman, 

For higher forces (/ > IQpN)^ the end-to-end distance L is close to Lq and the elastic restoring force is due to 
distortion of the internal structure of DNA. In this regime, the force extension curve can be approximated by 

/ ~ ksT-tDNA^^^r^, L>Lo (14) 

We will call 7 ~ {df /dL) Lo/ksT the "stretching modulus". 7dna is about 3GGnm~^ (Cluzel et al., 1996; Smith et 



al., 1996), i.e., almost four orders of magnitude larger than the corresponding value 3/Adna obtained from Eq. 13 



B. Bending and twisting of linker 

To calculate the stretch modulus of the {6, (j), b) model, each linker is modeled as a wormlike chain (WLC) of fixed 
length b (see Schlick, 1995 for a review of the WLC). We denote the geometrical configuration of the fcth linker 
{k = 1, 2, 3, ...) by rj, (s) with s being the arclength, < s < 6. The elastic energy of the wormlike linker is given by 
the sum of the bending and the torsional energies: 

Here k is the bending stiffness which is related to the persistence length Adna of (linker) DNA by k = ksTAuNA- 
Furthermore, 1/Rk (s) = |c?^rfe (s) /ds^j denotes the curvature of the A;th linker at the point s along its contour. The 
torsional angle of the linker is rjk (s) and the torsional stiffness is C. The positions (0) and (6) of the two termini 
of the kth linker coincide with the termini of the neighboring hnkers, i.e., Vk-i (b) = (0) and rfc (&) = r^+i (0). 
Furthermore, we assume that the entry-exit angles have the fixed value 7T — 9 independent of the bending and twisting 
of the linkers. This means that the unit tangents fulfill the condition cos {6) = t/. {b)-tk+i (0), with t/^ (s) = dr^ (s) /ds. 

(i) Enthalpic elasticity: We study first the stretching of the planar zig-zag pattern ((j) = tt, 9 arbitrary). The 
undeformed zig-zag fiber is depicted in Fig. 5(a). In order to give a more accurate description of the mechanical 
properties of the fiber, we assume that the nucleosome core particles are not located at the sites where two linkers 
come together but rather slightly displaced, forming a stem configuration as shown in Fig. 5(a). This is the geometry 
obtained from the electron cryomicrographs of Bednar et al. (1998) and it is the same geometry that was adopted in 
the computer simulations of fiber stretching (Katritch et al., 2000). In the following we denote the actual linker length 
by b in order to distinguish it from 5, the distance between neighboring nucleosomes. For symmetry reasons, there is 
no torque on the structure, so that dr/k/ds = 0. The stretching of the fiber is achieved by a bending of the linkers 
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with the entry-exit-angle 6 remaining constant, of. Fig. 5(b). This leads to a deformation where the tangent vectors 
tfc (0) and tfc (6) remain parallel but undergo lateral displacement. We assume a small deformation of the linker with 
displacement u (s) from the straight configuration small compared to b. (Since u (s) is the same for all the linkers, we 
drop the index k from here on.) 

From the minimization of Ek, Eq. we obtain the Euler-Lagrange equation d^u/ds^ = 0. The boundary 
conditions that must be obeyed by the solutions are u (0) = u' (0) = u' (b) = and u (Fj — d where d describes 
the displacement of the bead vertical to the original straight linker (we assume d <^h here and neglect terms of the 
order {d/b) ). It follows that the deformation profile is given by u (s) = —2ds^/r + Zds^ . The associated bending 
energy is E — Gnd^/b'^ (per linker). The deformation translates into an effective change in the deflection angle from 
9 to 9 — A9 where A9/2 = d/b - see Fig. 5(b). The energy of a fiber with N linkers as a function of A9 is thus given 
by£;=(3/2) (n/b) (AOfN. 

The change in 9 produces a change in the overall length of the fiber. We find from Eq. 0: 

L = bN cos i — - — ] ^ Lo + bNsm{9/2)A9/2, < 1 d < fe) (16) 

where Lq is the contour length of the unperturbed fiber, Eq. |l|. The energy can be rewritten in terms of the extension 
AL = L — Lq. The restoring force follows then from / = dE/dL: 



f^TT-TTTTT^^T-AL (17) 



12 K 
Nsin^ (61/2) P 

The associated stretching modulus (defined as in Eq. |l^ follows from Eqs. |^ and |l^: 



UAuNA cos(g/2) 

Jfiber [9) ~ . 2//)/o^ (1^) 

0^ sm [9/2) 

We next consider the deformation of fibers with crossed linkers ((/> <C 1, 9 large). When such a fiber is stretched, 
linkers will be twisted as well as bent. Interestingly, a fiber with a high bending stiffness (k oo) can still be stretched 
just by twisting of the linkers. When one applies a tension to such a fiber each linker is twisted and the rotational 
angle changes by Acj) from one linker to the next, i.e., drj/ds — A(/)/6. The twist is distributed homogeneously along 
the linker since d'^rj/ds'^ — which follows from minimization of Ek in Eq. The energy per linker is given 

hy E = (C/2) (A0^/6). The twist of the linkers changes the length of the fiber, and using L from Table I it is 
straightforward to calculate the force as a function of the relative extension: 

In the opposite limit of linkers with extremely high torsional stiffness (C ~* oo), the force-extension curve can be 
mapped onto the planar zig-zag case, described above, by replacing tt — by (/) ~ (/)Cot(0/2), the "effective" angle 
between two consecutive linkers as seen from the "side" of the fiber. (This follows from ip ~ 2Al/b where Al is the 
difference in the longitudinal position of bead i and i + \; Al = {(j)b/2) cot {9/2)) Using Eq. |l^ we find 

/ , ^^Vt^^ (20) 

b^N cos^ U/2] 

We can now define as before the two stretching moduli jtwist and ^bend using Eqs. 19 and ^ together with Eq. 
If we allow both twist and bend then the two "spring constants" act "in series": jj^^ber — 1 twist + Ibend- obtain 
for the stretching constant of the fiber ((/> ^ 1, large) 



la A\ 111 Ibend \ 6AdNA 

1 fiber [9, <P) = 76erid 1 + 



Scot^ (I) K 



(i) 



C 



(21) 



For large n (with n defined as n = tt/ {-k — 9)) the twisting contribution can be neglected since then ^bend/l^ 



twist 



(3/4) (7r/n)^ (k/C). For DNA C ~ ksT x 750A> k ( Klenin et al., 1989; Crothers et al., 1992) so that -fbend < Itwtst 
for 71 > 5. In this case one has "f fiber — Ibend and Eq. |2^ applies. 
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(a) Entropic elasticity: Just as for naked DNA, the entropic contribution to the elasticity dominates for weak forces 
{L Lq). The restoring force is again of the form 



with A f iter the persistence length of the fiber. This persistence length is calculated in Appendix C. For the case of 
the crossed-linker fiber we find 

A/^ber«^DNA^COt(0/2) (23) 

For values of and </> appropriate for the crossed-linker structure, Afi^er is somewhat less than Adna (of the order 
of 50nm). This surprising conclusion is related to the fact that a large amount of DNA material is stored in the fiber 
per unit length. The 30-nm fiber is thus indeed highly flexible. 

Our calculation of the stretching properties of the two-angle model predicts an important difference between the 
stretching behavior of DNA and that of the 30-nm fiber. The enthalpic stretching modulus of a chromatin fiber is of 
the order A-q^p^/W (see Eqs. |l^ and ^l]) which is about O.Bnm^^ for linkers with 40bp's (5 = 40 x 0.34nm ~ lAnm). 
This is only an order of magnitude larger than the entropic stretching modulus l/A/if,er — O.OSnm^^ . In other words, 
because of the low value of "/fiber and because of A fiber ~ ^dna there is no longer a very clear distinction between 
entropic and enthalpic behavior as it is observed for naked DNA. In conclusion, the 30-nm fiber shows soft elasticity 
under stretching due to bending and twisting of the linkers. 



C. Internucleosomal attraction 



The effect of attractive interaction between nucleosomes is to cause a compression of the 30-nm fiber. Phase behavior 
studies of linker- free nucleosome solutions, i.e., solutions of disconnected nuclesomes (Livolant and Leforestier, 2000, 
cf. also Fraden and Kamien, 2000) indicate that nucleosome core particles spontaneously form fiber-like columnar 
structures, presumably due to attractive nucleosome-nucleosome interaction. Attractive nucleosome interaction could 
be mediated for instance by the lysine- rich core histone tails (Luger et al., 1997), as mentioned above. 

It is important to distinguish these condensed fibers from the swollen solenoid-, zig-zag- and crossed-linker structures 
predicted by the {6, 0, 6)-model. The dominant energy of the condensed structures is the nucleosome attractive 
interaction, while the "swollen" structures are dominated by linker elasticity. In this section we will discuss the 
competition between swollen and condensed phases for a simple case. 

For simplicity, we model the fiber as a planar zig-zag structure with elastic linkers and assume in addition a short- 
range interaction between nucleosomes. This interaction, denoted by Uinter, is assumed to be a short range attraction, 
of strength —Umim that acts only when the nucleosomes are in close contact, i.e., at a distance a; « 2a of the order 
of the hardcore diameter. For a given nucleosome, say the ith, the closest nucleosomes in space are number i-\-2 and 
i — 2 as discussed in Section 2. We will disregard the interaction between other pairs. The elastic interaction Uei 
follows directly from Eq. ^ with N — 2: 

Uhendix) = —^——r^ix- Xof ^ —{x~ Xof (24) 
sm (0/2) o-' 2 

where Xq = 2&cos (0/2) denotes the distance between nucleosome i and i -I- 2 for straight linkers (cf. Eq. The total 
internucleosomal U (x) equals Uinter (x) + Uhend {x). 

Fig. 6(a) shows U {x) for different values of 9. We assume for simplicity that the interaction energy Uinter remains 
unchanged. Curve "1" in Fig. 6(a)) shows U (x) for a small value of 9 where the global minimum of U (x) is located 
a.t X — Xq denoted by " S" (swollen state) . Curve " 2" corresponds to an intermediate value of 9 at which the minima 
at "S" and "C" have the same value. For this value oi 9, 9 = 9c, the energy minimmn shifts from "S" to a new 
minimum, representing the condensed state " C" . The change in 9 produced a structural transition from a swollen 
state to a condensed state. Finally, curve " 3" depicts U (x) for a defiection angle 9 > 9c with the minimum at " C" . 
The critical angle for the "S" to "C" transition can be determined by comparing the bending energy at close contact, 
Ubend (2a), and the strength Umm of the short range attraction. Equating both leads to the following condition for 
9c-- 
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COS {0^/2) ^ ^^^5^ sin {0^/2) = J (25) 

In the swollen state the elastic properties are those discussed in the previous section. In the condensed state, the 
elastic properties are determined by the detailed form of the nucleosome interaction potential. 

If the condensed state has a lower free energy, i.e. ii > 0c, then an external stretching force / can induce 
a transition from the condensed to the swollen state. The transition point fcs follows from a "common-tangent" 
construction. The conditions are U' (xi) — U' {X2) — fcs and {U (X2) — U {xi)) / {x2 — xi) — fcs (cf. Fig. 6(a)). 
The first pair of conditions leads to xi = 2a, X2 = xq + fcs/ K. The last condition leads to 

fcs = V2KU~ -K{xo~ 2a) (26) 

The corresponding force-extension curve has a "coexistence plateau", cf. Fig. 6(b). If the imposed end-to-end 
distance is smaller than Lq (the contour length of the condensed fiber) then the restoring force is entropic. For 
Lo < L < Li the force rises sharply with increasing L. This "hard elasticity" is governed by the nucleosomal 
interaction potential Umter- Then at L = Li the coexistence plateau is reached. Between L ~ Li and L = L2 parts 
of the fiber are in the "S" state and parts are in the "C" state. For larger extensions, L > L2, the fiber shows soft 
elasticity due to the bending (and twisting ) of the linkers as discussed in the previous section. 



D. Stretching chromatin 

We now compare the results of the previous two sections with the force extension curves found in recent experiments 
(Cui and Bustamante, 2000). 

(j) Low ionic strength: We start with the force-extension profile measured at low ionic strength (5mM), cf. Fig. 2 
in Cui and Bustamante. As discussed above, at low ionic strength the chromatin fiber constitutes a swollen fiber with 
crossed linkers. The nucleosomes are far apart and we assume that there is no direct interaction between nucleosomes. 
The resulting force-extension profile is expected to show a crossover between an entropic elasticity (cf. Eq. ^2|) and 
a soft enthalpic elasticity with a stretching modulus given by Eq. |2l|: 



s'^^'^pL for L < Lo 

2knT^f..„ ^L-Lo) + ^ forL»Lo 



ipb4.N" " ^" 



rfibN 



with Lq ~ (^/2^ 



bN, cf. Table 1. 



Cui and Bustamante estimate the number of nucleosomes in their fibers to be iV « 280. From the formula for 
L in Table I, we would estimate the length of the fiber to be Lo « 1.0/j,to using the values — 95° (Bednar et al., 
1998), 4) = 36° (Widom, 1992), and b = AObp = 40 x 0.34nm 1.4 x IQ-^m. The linker length is estimated from the 
nucleosome repeat length of roughly 210 bp's (cf. Table 7-lF in van Holde, 1989) minus roughly 170 bp's that are 
associated with the core and linker histones (cf. page 268 in van Holde, 1989). Using the moduli for DNA (Hagerman, 
1988), K = ksT X bQnm = 2 x lO-^^pNm^, C = ksT x 75nm = 3 x lO-^^pNm? and Ip = 30nm for the persistence 



length of the fiber (cf. Appendix C) we find from Eq. 27 the following force-extension relation (force inpN, extension 
in ^m): 

r 0.35 X L for L < 1.1 

- \ 1.2 X (L - 1.0) -f 0.4 for L > 1.1 

The agreement with the experimental curve at low ionic strength (5mM NaCl) is reasonable (cf. Fig. 2(a) and (b) 
in Cui and Bustamante, 2000). More explicitly, for forces up to 5pN and extensions up to « 2/im there are two 
distinctive regimes: For small extensions, L < 1/im, the force increases only slightly with tension, namely roughly as 
/ « 0.5 X L. Then for L > l/im the measured force increases much faster and shows the following linear dependence: 
/ w 7 X L (Fig. 2(a)) or / w 5 x L (Fig. 2(b)). The different slopes in this regime are a result of a slight hysteresis: 
the relaxation curve has a smaller slope after the fiber has been stretched to an end-to-end distance 2.5fj,m (Fig. 
2(b)) than for the case of a much smaller stretching cycle (up to 1.8/im, Fig. 2(a) in Cui and Bustamante, 2000); the 
hysteresis disappears for a smaller rate of extension or contraction, and might be a result of nucleosome-nuclcosome 
interaction or of modifications of the fiber close to the entry-exit point of the linkers at higher tension. We also 
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mention that for forces beyond « bpN the (relaxation) curve shows an increasing sfope, probably due to nonlinear 
effects not accounted for in the current study^. 

The calculated forces are smaller than the measured ones (roughly by a factor of 4), for several reasons. First, the 
(mean) values of 0, (j) ^^nd h (and thus TV) are only roughly known. Secondly, the value of we used (95°) is not 
large enough compared to tt for the above given theoretical formulas to hold accurately. However, as a check of our 
analytical approximations, we compared our results with the computer simulations by Katritch et al. (Katritch et 
al., 2000) where 6*, (/) and h are variable. This comparison is given in Appendix D, where we show that there is good 
agreement indicating that our analytical approximations were in fact reasonable. 

{a) High ionic strength: For 40mM NaCl or higher ionic strength the chromatin fiber is much denser and nuclcosomes 
approach each other closely. Attractive short-range forces and the increase of associated with higher ionic strength 
should favor the condensed phase. A plateau indeed appears at bpN in the force-extension plot (cf. Fig. 4 in Cui and 
Bustamante, 2000). From the extent of the plateau, 0.6/j,to, its height, 5piV, and the number of nucleosomes in the 
stretched fiber, w 280, it was estimated that there is an attractive interaction energy of roughly SfcT per nucleosome 
(Cui and Bustamante, 2000). 

We now can use Eq. ^ to estimate independently the strength of the nucleosomal attraction from the value of the 
critical force alone. We find: 

_ {fcs + K{x^~2a)f 

^ min — 2_/^ V / 

If we neglect the second term in the bracket, we find Umin ~ /^/ (2i^) ~ 6fcT (assuming 6 — 140°), close to the value 
3fcT estimated directly from the force-extension diagram (Cui and Bustamante, 2000) and also in accordance with 
the computer simulation of Katritch et al. who obtained an internucleosomal short-range attraction of order 2kT 
(Katritch et al, 2000). 

Using Umin = ikT we can estimate the critical value 9 — 9c s^t which the condensed and the swollen chromatin 
fiber should coexist. We find numerically from Eq. |2^ that 9c ~ 100° (using a = 5nm and b — 14nm). This value 
is lower than the one that can be inferred from experiments. At ISmAf NaCl {9 « 135°) the fiber appears to be 
decondensed, as indicated by stretching experiments and from electron cryomicrographs. This fact as well as the 
appearance of a plateau in the force extension curve at 40toM salt (where 9 w 140°) indicates that one should expect 
135° 9c ^ 140° (cf. Bednar et al., 1998). It should be recalled, however, that our model for the attractive interaction 
is highly oversimplified. 



V. CONCLUSION 



The present analytical study of the {9, cj), b) model first of all shows that this model can account for the measured 
force-extension curve of the 30-nm fiber in the low-salt regime with, in effect, no fitting parameters (since 9, cf), and 
b can be estimated experimentally and since the elastic moduli characterizing naked DNA are known). Since the 
{9,(j),b) model also accounts for the observed low-salt structure of the 30-nm fiber ("crossed linkers"), there seems to 
be good evidence that this model is at least the proper description in the low-salt regime. 

We have been able to compute the structural and elastic properties over a wide range of {9, 0)-values. We suggest 
that the native chromatin fiber might be a particular realization of this rich array of structures, namely the one that 
simultaneously maximizes compaction and accessibility, consistent with the restriction of excluded volume between 
nucleosomes. 

Confirmation that a certain optimization principle is in fact operative for biomolecules is usually a difficult issue. 
We already saw that, at best, the principle is incomplete since the linker-length b evidently is not determined by 
the conditions of maximum compaction and accessibility. One possibility may be to explore the fine-structure of the 
dotted curve in Fig. 4, the lower bound of as a function of 9 . This is expected to have an "irregular" shape due 



Our calculation is based on the assumption of small deformations. For the zig-zag case this requires d ^ b, i.e., A9 <^ 1, 
cf. Fig. 5(b). Using Eqs. [| and this condition translates into the requirement that the tension / is smaller than 
6k/ (P sin (61/2)). For fib ers with internal linkers the condition is A(f> 1, leading to / <g 6k/ cos^ ^0/2^^ . Thus in both 

cases a good estimate for the range of forces where the linear approximation holds is given by / < 6k/P. For the chromatin 
fiber under consideration we find 6k/&^ ~ 6pN. 
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to commensurate-incommensurate effects and it may be possible to associate a discrete geometrical structure (e.g. a 
particular index n for the polygonal star projection) with maximum compaction and accessibility. Such a study would 
require, however, a better description of the structure of individual nucleosomes and extensive numerical work. 

How confident can we be that the {9, (j), b) model is appropriate as well in the biologically relevant regime of physi- 
ological salt concentrations? We had to include a weak attractive nucleosome interaction to explain the coexistence 
in the force-extension curve. If the fitted value for the attractive potential [Umin) is used in Eq. ^ we obtain a 
reasonable estimate for the critical angle Be for the "S" to "C" transition (but with a significant error). 

A completely different approach would be that the high-salt regime is controlled not by a balance between soft 
elasticity and weak attraction but completely by nucleosome- nucleosome attraction forces (plus short-range repulsion) . 
As shown by the work of Livolant and Leforestier (2000), nucleosome attraction indeed can produce discoidal fiber 
structures (formed by linker-free core particles) all by itself. If the interaction energy is strong enough, then the 
linkers would be strongly bent in the condensed state. The {9, (p, b) model would not be a valid description anymore. 
The effect of tension could be to produce a sequence of different condensed structures. Only at high tension when 
the internucleosomal contacts are broken one recovers the soft-elasticity regime, described well by the {9, 0, b) model. 
Which of the two approaches is valid is an issue that must be determined experimentally. 

Interesting questions for "chromatin physics" in the future may focus on dynamical issues. Suppose that 9 is locally 
increased, e.g. by acetylation of core histone tails, how long does it take for the accessibility to increase sufficiently. 
How important is nucleosome mobility (Schiessel et al., 2000) and nucleosome "evaporation" (Marko and Siggia, 1997) 
for the swelling dynamics of chromatin? 
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APPENDIX A: THE MASTER SOLENOID 



For any given set of angles (0, 0) there is a solenoid so that the successive monomers of the fiber structure lie 
successively on this helical path. (There are actually many such solutions, but we are interested in the one with the 
largest pitch angle 7.) We parametrize the solenoid as follows 

(i? cos (as/R) \ 
Rsm{as/R) (30) 

R denotes the radius of the solenoid and a is related to the pitch 7 by 

a = cot 7 (31) 

(as follows from r (0) = (0, a, 1)). 

Assume now an infinite fiber of monomers with a given pair of angles (0,0). The monomers are located at the 
positions Ro,R±i,R±2, •■• The axis of the fiber coincides with the Z-axis. Assume further that we choose the values 
R and a so that the solenoid curve goes through all monomers. Put the monomer labeled i = at s = so that 
Ro = {R, 0, 0); the subsequent monomer, i = 1, is at a position Ri given by Eq. ^ with s ^ sq. The next monomer 
is located at R2 = r {2sq). Finally, the position of monomer i = — 1 is given by R_i = r (— so)- 

Now let us calculate the bond vectors between these monomers. Monomer i = 1 is connected to monomer i = via 

R cos {aso /R) — R 
ro = Ri — Ro = I i?sin (aso/i?) 

So 

The separation vector between monomer i = 2 and i = 1 is given by 

(R (cos (2aso/i?) — cos (asa/R)) \ 
R (sin (2aso / R) — sin (asQ /R)) 
So J 

and that between monomer i = and i = — 1 by 

/ R- Rcos{aso/R) \ 
r2 = Ro — R-i = i?sin(aso/-R) 

\so J 

So follows from the condition of fixed linker length, i.e., |ro| = b. This leads to the relation 

= 2R^ (1 - cos (aso/i?)) + si (32) 
We determine from cos = ro • r2 / | | , which leads to 

2R^ cos [aso/R] {I- COS. {aso/R)) + si 

cosO = TTT^TT. rrrr^ 5 (•j-j) 

2R^ (I - cos (aso I R)) + si ^ ' 

Finally, is the angle between normal vectors of the planes that are defined by monomers and 1, i.e. cos ni • n2. 
We obtain ni and w.^ from ni = A/ |A| and n2 = B/ |B| where A = ro x ri and B = r2 x ro. After some algebra we 
arrive at 

cos (aso/i?) + i?^ sin^ (aso/i?) 

cos0= T. 9- ■ — ^ (34) 

sg + i?2 sin^ (aso/i?) 

Equations p3, and p3 relate ol (or 7), R and so of the spiral to 0, B and b. 



APPENDIX B: RANDOMNESS IN THE 0-DISTRIBUTION 
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Up to now we have assumed that the values of the angles 9 and (j) are constant throughout the fiber. The resulting 
"ground state" configuration (unbent and untwisted linkers) is a fiber whose axis is perfectly straight. The assumption 
that the linker entry-exit angle 9 is constant is based on the fact that it is a local property of the nucleosome core 
particle, and as long as the biochemical conditions are homogeneous throughout the fiber this should be a reasonable 
assumption. It is known, however, that the rotational positioning is not perfect, as can be seen from the experimentally 
determined distribution of the linker length in chromatin (Widom, 1992). Even though a preferred rotational setting 
can be deduced, the width of the distribution of linker lengths will be reflected in the width of the distribution of 
the angle 4>. If the rotational setting of the nucleosomes were completely random, then the chromatin configurations 
would correspond to particular configurations of the freely rotating chain (if we neglect excluded volume effects) 
(Doi and Edwards, 1986). These configurations, in turn, are those of a Gaussian chain with a persistence length 
Ip = 6(1 + COS0) / (1 — COS0) (the Kuhn statistical length as defined in Doi and Edwards, 1986). Note that Ip 
increases when 9 decreases, a mechanism similar to the accordion-like unfolding of the zig-zag structure or the 
untwisting of the fiber with crossed linkers discussed above. In the following we will assume small variations of the 
rotational setting around some mean value 0. We consider the three cases: the solenoid, the fiber with crossed linkers, 
and (twisted) zig-zag structures. 

(?) Solenoids (iji <C 1, <C 1): We start with the solenoidal fiber with ^ ^? <C 1. Then the pitch angle is small (cf. 
Eq. ^) and each loop of the solenoid resembles nearly a circle. The small variations in (j) will add up to an effective 
deviation from the original orientation of the fiber per turn of the helix. If one has n monomers per turn it can 
be shown that (A^^) = ntr^ with the width of the (/)-distribution. A^ is Gaussian with a width (jq — \pnu^. With 
each turn the middle axis of the solenoid proceeds by a length d where d is given by Eq. ^j. We can interpret the 
middle axis of the solenoid as a new effective chain with bond length d, and calculate the average of the scalar product 
of an arbitrary pair of successive bond vectors and a^-i-i of this new effective chain: 



27ra^ 



a,a,+i) = / dACcos (AC) exp \ ^ d' \ 1 - ^ \ , (35) 



the approximation holding for <^ I. It follows that the end-to-end distance of the chain, assuming that the solenoid 
has AI turns (corresponding to a fiber of = 2ttM/9 nucleosomes), is 

MM , ,9 It r 

) = E E (^^a™) ^ (36) 

n— 1 m— 1 C 

The persistence length Ip of the fiber follows from Ip = {L'^) / {dM)\ 

IP . ^6 (37) 

where we have made use of the relations d ~ (27r0/6'2) h (cf. Eq. ^ and — 2'Ka'^/9. These results must be modified 
for solenoids with larger pitch angle 7 (0 <C ^ 1), where - since only the component A0cos7 of a variation A(/) in 
the rotational angle (j) leads to a change in the direction of the fiber (the component A0sin7 leads to a twist) - one 

has to replace by cr^ cos 7. The resulting persistence length is given by /p ~ 46/ {a'^ cos^ 7^ ~ AbcfP' / (cr^^^) ■ 

(ii) Fiber with crossed linkers {(p <C 1, tt ^ 9 <C 1): This case can be calculated analogously. The number of 
monomers per "turn" is given by n = tt/ (tt — 9) (see above) so that — y/ir/ (tt — 9)a^. Furthermore, the bond 
length of the new effective chain is d ~ ii(f}h/A (cf. Eq. From Eq. |3^ it follows that a fiber of = nM monomers 
has the mean-squared end-to-end-distance 

(i^) ^ llS^b^M (38) 



and a persistence length 



I, ^ il^, (39) 



Now consider typical values for chromatin: 36°, h — 20nm, and 9 « 145° at 80mAf, 9 « 135° at 15mM and 
9 « 95° at 5mM (Bednar et al., 1998). Assume that the histones are located at equidistant positions but with small 
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variations, typically ±l&p, i.e., « 36°; then we find Ip « 20nm at 80mM, Ip « 25nm at 15mM and Ip « Alnm 
at 5mM. 

{in) Twisted zig-zag fiber. Finally, we consider zig-zag structures, first the case where the angle (p fluctuates around 
the mean value tt (planar zig-zag). (With no fluctuations in the zig-zag structures simply represent a perfectly 
flat ribbon.) Assume first that one bond is slightly rotated by A(p <C 1. As a consequence the ribbon is deflected 
by an angle A^i ~ sin (6/2) Aip; furthermore the orientation of the plane defined by the ribbon rotates by an angle 
~ cos {6 /2) Acp. A long ribbon-like zig-zag structure with small fluctuations of the 0-angle shows individual 
configurations typical of a polymer with an anisotropic bending rigidity (Nyrkova et al., 1996); such a polymer has 
a plane of main flexibility, being highly rigid in the direction perpendicular to this plane. The anisotropy leads to 
two persistence lengths; an in-plane persistence length li which is associated with the deflection of the ribbon within 
the plane of main flexibility; and an out-of-plane persistence length I2, the typical polymer length that is needed to 
"forget" the orientation of the plane of main flexibility, li follows from the number of monomers ni that is needed on 
average to forget the original orientation of the axis of the ribbon, A^fni = 4 (the numerical value is chosen so that 
li is compatible with the definition of the Kuhn statistical length). Thus 



4 cos (61/2) b 
sin^ (61/2) 



h = bcos (6/2) m ~ "7 \ -3 (40) 



Similarly, I2 follows from (27r)^ = AC|n2: 



We consider the two limiting cases, (i) 9 = 0: here h = 00 and I2 — ^ib/a"^. The configuration of the chain is that 
of a straight line. Variations in A(p do not affect the positions of the monomers, (ii) 9 n: By reaching this limit 
the chain collapses into a configuration where it just goes back and forth between two monomer positions. Indeed, 
we find from the above equations that Zi — > and I2 c>o. 

For a twisted zig-zag structure with (jj = n — S with 6 1 there is an inherent orientational persistence length that 
follows from the twist of the fiber. This leads to a length I2 — bcos (9/2) n2 where fi2 = 2tt/6 denotes the number of 
monomers per turn. Apparently the inherent twist competes with the randomly introduced one, and the out-of-plane 
persistence length I2 (cf . Eq. ^l|) has to be replaced by I2 if h '^h- The role of variations in the linker length in the 
case of a twisted zig-zag structure was simulated by Woodcock et al. (cf. Fig. 3 in Woodcock et al., 1993). They 
chose the case 9 = 120° and S ^ 360°/13 ~ 0.48. Using Eqs. |^ and |l| we find h ~ 2.7b/al, I2 ~ 8.06/cr^ and 
I2 = 6.56. It follows from our formulas that the persistence lengths li and I2 decay rapidly with a^, a trend that can 
also be seen in the displayed configuration in Fig. 3 of Woodcock et al. If we choose, for instance, cr^ = 1/2 we find 
Ip ~ 116, a persistence similar to that of the fiber displayed in their Fig. 3b. If we double a^, i.e., cr^ = 1, we find 
li ~ 36 so that there is no longer a well-defined fiber; a similarly disordered fiber is displayed in their Fig. 3d. A 
closer comparison between our theoretical results and the disordered fibers shown in Woodcock et al. is not possible 
since in the case of the " simulated" fibers a discontinuous distribution of the values of (p was chosen, thereby varying 
the number of base pairs per linker. 



APPENDIX C: PERSISTENCE LENGTHS 

We calculate here the effect of linker flexibility on the persistence length of the two-angle fiber. We first calculate 
the zig-zag-structure where one has two different persistence lengths, the persistence length /p"-* for bending within 
the plane of the fiber, and the length Ip^*^ for bending out of the plane. 

(j) bending in the plane of the fiber: Assume that the ribbon is bent within its plane with a large radius R of 

curvature so that R ^ b. The linkers are bent but not twisted in this case. Up to corrections of order ib/R) the 
shape of each linker (i.e. its deviation from a straight line) is given by u (x) = —ex^ jb -\- ex. This function fulfills 
the appropriate boundary conditions w (0) = u (6) = and u' (0) = —v! (6) = e. This leads to the following bending 
energy per linker: E = (K6/2i?^) cos^ {d/2). In the longitudinal direction of the fiber this corresponds to the bending 
of a piece of the length 6 cos {9/2). Thus 

Afl,=A^^p,cos{9l2) (42) 
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{a) bending perpendicular to the plane of the fiber: This bending is accomplished by a combination of twist and 
bending of the hnkers. Consider the two cases separately. If there is only twist allowed (k oo), then each linker 
has to be twisted by an angle 6cot (9/2) /R which leads to a twisting energy E = C cot^ {0/2) b/2B? and then in turn 
to the persistence length = {C/ksT) cos {9/2) /sin^ (^/2). Now consider the case without twisting {C oo) 

but with bending of the linker only. If one bends a linker out of the plane of the fiber with a radius of curvature 
R, it can be shown that as a result the zig-zag is deflected by an angle fecos {9/2) / R. If each linker is bent in such 
a way, the zig-zag fiber is bent out of its plane with an overall curvature of 1/R. The bending energy per linker is 
E/ksT — Adna&/2-R^, leading to a persistence length A^^^J — Adna/ cos (0/2). By putting the two deformation 
modes "in series" we find the overall persistence length for bending the zig-zag out of the plane: 



^ , -I ^ _ ^DNA 1 

S^ber " \^l^t^^st' + ' I ^bend ) ^ (^/2) 1 + % tau^ {9/2) ^ ' 



9/2)': 

(out) 



For small angles of 0, the bending contribution dominates and Af^^^J^ ^dna for ^ (naked DNA). On the 
other hand, a very dense zig-zag with a value of 9 close to tt is bent by the twisting of the linkers, leading to a very 
short persistence length — {C /ksT) cos {9/2). Interestingly, for DNA where k C one finds from Eq. ^ 

^/°i6er ~ ^DNA COS {9/2) Over the whole range of 0-values. Thus, in this case A'^^^J^^ « ^fib2-- 

We turn now to fibers with crossed linkers {(j) small, 9 large). If we bend such a fiber within a given plane, then an 
inhomogeneous deformation pattern result where some of the linkers are oriented (nearly) parallel to the fiber while 
others are perpendicular. The first class of linkers will be bent, the second will be mostly twisted. The effective angle 
is now TT — 4>cot {9/2) instead of 9 (this follows from L in Table I with N = 2). Since C ~ k the contribution to the 
elastic energy is approximately the same for all the linkers, leading to a persistence length 

Afaer - ^DNA COS ~ cot {9/2)^ ~ Aqna ^ cot {9/2) (44) 

Using the 6'-values given by Bednar et al., k, 36° (and ^dna ~ SOnm) we find Afn,er ~ 14nm for 9 = 95° (the value 
at 5mM monovalent salt), Afn,er ~ 6nm for 9 = 135° (15mM) and Afiber ~ 5nTO for 9 — 145° (80mM). These values 
are smaller than the diameter of the fiber so it is reasonable to assume A fiber ~ 30nm for the persistence length of 
the 30-nm fiber (roughly independent of the salt concentration) . 

APPENDIX D: COMPARISON WITH COMPUTER SIMULATIONS 

Katritch et al. performed Monte-Carlo simulations of the two-angle model with flexible linkers (Katritch et al., 
2000). Here we compare their results with our theoretical predictions. We first base our analysis on the results for 
the zig-zag case, Eq. |l^, for reasons given below. Including the entropic contribution we find the following force law: 



f ^ J /Jofw2) ^ fo'^ ^ « COS {9/2) 

^ ihT^ ^ - cos {9/2)] + 3^ for x » cos (0/2) ^ 



To allow a better comparison with the diagrams in Fig. 3 of Katritch et al. we present in Eq. ^ the force as a function 
of the relative extension x ^ L/ (bN) . While most of the data of Katritch et al. were obtained assuming a (quenched) 
set of random values of the rotational setting of the nucleosomes, we believe that the qualitative dependence of / 
on 9 and b should be unaffected by this assumption. Figure 3(a) in Katritch et al. shows the dependence of the 
force-extension profile on the entry-exit angle W ^ tt — 9. It can be seen that the initial slope (entropy regime) 
decreases with W (increases with 9) in accordance with Eq. The behavior at larger forces shows the opposite 
dependence, as is also predicted by Eq. Finally, the crossover is shifted with increasing W to larger values; this 
again is in accordance with Eq. Figure 3(b) in Katritch et al. depicts the dependence on the linker length b; here 
the data show no indication of a dependence of the initial slope on b, and similarly for the value of the crossover, all 



in accordance with Eq. 45, Then for the second regime they find an increasing slope with decreasing linker length, 
also in accord the above theory. 

Fig. 3(c) of Katritch et al. shows a comparison between random rotational settings of the nucleosomes and non- 
random settings. The zig-zag case, cj) = tt, has a slightly smaller initial slope than the random case and a slightly 
larger slope in the second linear regime, but follows always quite closely the case of random settings. This justifies 
the above comparison between Eq. |4^ — based on the zig-zag case — and the computer simulations using a random 
setting. Figure 3(c) of Katritch et al. allows now also a direct quantitative comparison. We predict from Eq. E^, for 
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the values 9 = 2.27 and b = 40bp (used in Katritch et al.), that / ~ 0.95x for x < 0.42 and / ~ 14.9 {x - 0.42) + 0.4 for 
X > 0.42 (force in pN), in good agreement with their data points. (For x > 0.6 the datapoints indicate an increasing 
slope, a result of nonlinear effects not taken into account in our theory.) 

Finally, Katritch et al. also provide data points for the case 4> = 0.35 and 8 = 2.27. We find from Eq. |2^ that one 
has / ~ 4.9xforx < 0.1 and / « 8.6 (x - 0.1) + 0.5 for x > 0.1. (Here x ee L /67V, and xq = Lq/IN = (0/2) cot (6'/2) ~ 
0.1.) This result overestimates the entropic contribution as can be seen by comparison with the Tw = 20°-curve in 
Fig. 3(c) of Katritch et al. (data points are missing for x < 0.2 but the force at x = 0.2 is smaller than IpN). The 
reason is probably the underestimation of the persistence length of the fiber Ip, assumed here to be of the order of 
the fiber thickness. The real value could be larger since this set of angles corresponds to an extremely dense fiber 
where excluded volume effects become rather important. The simulation data for this set of angles is not found to be 
in good agreement with the experimental force-extension characteristics. This might be attributed to the small value 
of the entry-exit angle (130° instead of 85° as suggested by the electron cryomicrographs, Bednar et al., 1998); if we 
use 9 = 95° and = 36° (as above) we find / ~ 1.4x for x < 0.3 and / « 4.7 (x - 0.3) + 0.2 for x > 0.3 which is 
closer to the experimental curve. 
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Table I 
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crossed linkers 
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twisted zig-zag 
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1 l + ^^(tan^(0/21/.l-l)/2 


TTtjiU-b'-'' 
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Ta6le caption: 

Table I: Geometrical properties of the two-angle fiber for the three limiting cases: solenoid, fiber with 
crossed linkers and twisted zig-zag fiber. Displayed are the fiber radius R, the length L of a, {N + l)-mer, 
the line density A = N/L, and the 3D density p = X/ttR"^. 
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Figure captions: 

Figure 1: Schematic representation of the nucleosome. Eight core histones aggregate into the histone 
octamer that acts as a cyhndrical spool around which the DNA is wound in l-and-3/4 turns. The hnker 
histone is also depicted that acts at the entry-exit point of the DNA. The entry-exit angle tt — ^ of the 
linker DNA is one of the angles defining the two-angle model. 

Figure 2: The two competing models for the 30-fiber: (a) the solenoid model and (b) the crossed linker 
model (see text). 

Figure 3: Fraction of a two-angle fiber containing four nucleosomes (it is a part of structure "11" in Fig. 

4). The two angles arc depicted, the deflection angle 9 and the rotational angle (f>, together with the 
"nucleosome diameter" 2a and the "linker length" b. All four are considered to be constant throughout 
the fiber. The arrows denote the nucleosomal axes, cf. Fig. 1. 

Figure 4: Diagram of geometrical states of the two-angle model. Shown are examples of different configura- 
tions and their location in the (^, ^)-space (arrows). The dashed and dotted curves depict the boundaries 
to the {0, (/))-values that are forbidden due to excluded volume interaction, one regime (large 6'-values) 
due to "short-range" interaction, the other (small values) due to "long-range" interaction (see text for 
details). 

Figure 5: Stretching of a zig-zag chain, (a) The unperturbed chain {F = 0) has a total length Lq and 
straight linkers, (b) The same fiber under tension F > 0. The fiber is stretched to an end-to-end distance 
L > Lq hy bending of the linkers. The linkers are bent in such a way that the entry-exit angles at the 

individual nucleosomes remain unchanged. 

Figure 6; (a) Internucleosomal interaction potential U between nucleosome i and i -I- 2 as a function of 
distance x. In addition to the elastic contribution there is an short range attraction for nucleosome at 
close contact, x = 2a. The different curves correspond to different values of the angle 6. Curve "1" has 
the global minimum at large x (swollen state "S") whereas curve "3" has the minimum for nucleosomes 
in close contact (condensed state "C"). Curve "2" corresponds to the transition point. Also depicted is 
the common tangent for curve "2". Its slope corresponds to the critical stretching force fcs at which 
nucleosomes are transferred from the condensed to the stretched state, (b) Force-extension curve of a 
condensed fiber (say, the fiber with the interaction potential "3"). For extensions L with L\ < L < L2 
one finds a coexistence plateau with the restoring force fcs (see text for details). 
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